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5} : ABSTRACT 

Studying the average properties of active galactic nuclei (AGN) host stellar 
populations is an important step in understanding the role of AGN in galaxy 
evolution and the processes which trigger and fuel AGN activity. Here we calcu- 
late model spectral energy distributions (SEDs) that include emission from the 
AGN, the host galaxy stellar population, and dust enshrouded star formation. 
Using the framework of cosmic X-ray background population synthesis modeling, 
the model AGN hosts are constrained using optical (B band) and near infrared 
(J band, 3.6 /im, 5.7 /zm, 8.0 /im, and 24 /im) luminosity functions and number 
counts. It is found that at z < 1, type 1 and type 2 AGN hosts have similar 
stellar populations, in agreement with the orientation based unified model and 
indicative of secular evolution. At z > 1, type 2 AGN hosts are intrinsically 
^ ! different from type 1 AGN hosts, suggesting that the simple orientation based 

unified model does not hold at z > 1. Also, it is found that if Compton thick 
(CT) AGN evolve like less obscured type 2 AGN, then, on average, CT AGN 
hosts are similar to type 2 AGN hosts; however, if CT obscuration is connected 
to an evolutionary stage of black hole growth, then CT AGN hosts will also be 
in specific evolutionary stages. Multi-wavelength selection criteria of CT AGN 
are discussed. 

Subject headings: galaxies: active - - galaxies: quasars: general - - galaxies: 
Seyfert — galaxies: stellar content — X-rays: diffuse background 



1. Introduction 

I t is known that all massive g alaxies have a central supermassive black hole (SMBH) 



(e.g., 



Kormendy fc Richstondll995[ ). Yet it is unknown why only a small fraction of these 
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SMBHs are actively accreting as active galactic nuclei (A GN). Several theories attempt to ex- 



plain the fueling mechanism of AGN: major mergers (e.g.. lSanders et al.l ll988; Hopki ns et al. 



2003 



2006) , minor mergers and gravitational instabilities within the host galaxy (e.g . , I Crenshaw et al 



bursts 
bursts 



e.g., 



Davies et al. 



Chen et al. 



2007; 



2009 



with the nuclear region (IMcKernan et al. 



Kormendv &: Kennicuttl 12004 ; Pierce et al.l l2007t IShapiro et ajj [2008 ), nuclear star- 



Ballantvnel 20081). supe rnova explosions within nuclear star- 



Kumar fc 



ohnsonl 120101 ) , and collisions of warm halo clouds 



2010b|). Given the large range of observed AGN 



properties, it is likely that different fueling mechanisms come into play for different AGN 
populations. Current observations suggest that powerful quasars are triggered by major 
mergers but moderate luminosity AGN are more likely to be triggered by stochastic fuel- 



ing in cidents ( IBallantyne et al.ll2006ai lHasingerll2008l ; iHopkins fc Hernquistll2009l ; iLutz et al. 
2010J). However, it is not yet clear which fueling mechanisms are dominant in which portions 
of the AGN population. 



The unified model of AGN (jAntonuccil 119931 ) explains the different observed levels of 
AGN obscuration as a simple geometric effect. This model assumes that the central engines 
of all AGN are identical and the level of obscuration is dependent on the line of sight between 
the observer and the central engine. Thus type 2 AGN are obscured by column densities 
Nh > 10 22 cm -2 , because the observer is looking through the dusty torus, in contrast to 
type 1 AGN which are unobscured (Nh < 10 22 cm -2 ) because the observer is looking down 
the throat of the dusty torus. Since the central engines are identical, it is expected that, 
on average, the host galaxies of various spectral types of AGN will be similar. However, 
different black hole fueling mechanisms are likely to lead to different relationships between 
the various spectral types of AGN. In the major merger AGN fueling paradigm, the AGN is 
triggered while deeply embedded in gas and dust and thus when the AGN first turns on, it is 
highly obscured. Eventually the radi ation pressure pushes away the rema i ning gas and dust 



revea ling an unobscured quasar (e.g., [Page et al.ll2004l ; IHopkins et al.ll2006l ; iRigopoulou et al. 



20091 ). In the major merger paradigm, type 1 and type 2 AGN have host galaxies which are 
in different evolutionary stages. Conversely, in the nuclear starburst fueling parad igm, type 1 
and t ype 2 AGN hosts can be similar since on average the unified model can hold (IBallantyne 
20081 ). This shows that tests of the unified model can be used to explore which processes 
are viable options for AGN fueling and in which sections of the AGN population different 
mechanisms are relevant. 

Classifying the host galaxies of Compton thick (CT) AGN, AGN with N H > 10 24 
cm -2 , is of special interest. It is expected that a large fraction of AGN with CT lev- 
els of obscuration are in a young evolutionary stage characterized by rapid black hole and 



host bulge growth dSanders et al 



1988 



Fabian 



1999 



Hopkins et al.ll2006l : iFabian et al.ll2009 



Draper &: Ballantynei I201CH : iTreister et al.ll2010l ). Thus, understanding the host properties 



-3- 



of CT AGN may offer special insight into the AGN triggering process and the formation of 
massive galaxies. Large samples of CT AGN are difficult to identify because, b y definition, 



CT AGN suffer from extreme levels of obscuration (e.g.. iGhisellini et al.lll994j ). Therefore, 
previous studies of AGN host galaxies have not been able to study the stellar populations of 
CT AGN. 

Galaxy optical colors are often described in terms of the red sequence and the blue cloud, 
where objects located on the red sequence are characterized by massive older stellar popula- 
tions and objects on the blue cloud are less massive and have young stellar populations and 
current star formation. Between the red sequence and the blue cloud is a lightly populated 
region referred to as the green valley. Some studies find that AGN feedback is connect ed to 
the shut down of host star formation (e.g., iMcKernan et al.l 1201 Oat iBluck et al.ll201l|) and 



thus certain types of AGN hosts preferentially r eside in the green valley ( e.g., iHickox et al. 
20091 ; iGriffiths fc Sternl boiol : IWold et al.l boiol : ISmolcic fc Riechersl J201lh . Thus AGN ac- 
tivity is potentially a stage in galaxy evolution where feedback from the AGN could shut 
down star formation, causing the host galaxy to age from the blue cloud ; across the green 
valley, and onto the red sequence (IFaber et al.ll2007l ; ISchawinski et al.ll2007l ). However, other 
studies find that the colors of an AGN and its host are more closely related to the amount of 
available obscuring material rather than the evolutionary stage of the host stellar popul ation 



(IBrusa et al. 



2009 



Georgakakis et al.ll2009l ; ICardamone et al.ll2010l ; iRovilos et al.ll201ll ). In- 



Li 



deed, ICardamone et al.l (120101 ) show that AGN are intrinsically bimodal in color and that 
many AGN which appear to be on the red sequence are actually located in the blue cloud 
once dust extinction is taken into account. Thus further investigation is necessary to es- 
tablish the nature of AGN host galaxy stellar populations and the role of AGN feedback in 
regulating star formation. This point is especially salient as several possible AGN fueling 
mechanisms include either causal or concurrent star formation. Some studies even show that 



AGN obscuration may in part be du e to nuclear starburst disks (e.g.. lThompson et al.ll2005 



Davies et al.ll2007t lBallantyndl2008l ) . In order to understand how host galaxy processes affect 
the central SMBH and the role of the central SMBH in host galaxy evolution, the nature of 
stellar populations of AGN hosts must be well understood. 

As AGN hosts exhibit size-able object-to-object variability, it is difficult to elucidate 
trends in the AGN population by fitting individual objects. Thus, in order to study larger 
trends in AGN-host galaxy interactions, the average stellar properties of a large ensemble of 



AGN hosts is investigated here. To this end the stellar synthesis models of lBruzual fc Chariot 



(120031) are used to explore the stellar populations of AGN hosts. As in previous stud- 
ies (jBallantyne et al.ll2006bt Draper fe Ballantvnell2011l). t he AGN SEDs are calculated us- 
ing the photoionization code Cloudy ( IFerland et al.lll998l ). The cosmic X-ray background 
(CXRB) synthesis modeling framework is used to characterize the AGN population. The 
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model SEDs are then used to move the CXRB framework into other wavelength regions. 
This allows the space density and evolution of AGN host galaxies to be determined by the 
most comprehensive census of AGN activity. By comparing these models against various 
observations in the optical through mid-infrared (mid-IR) spectral ranges, constraints are 
placed on the average stellar populations of type 1, Compton thin type 2 (here referred to 
simply as type 2), and CT AGN at various redshifts. In Section [2] the AGN spectral model 
and stellar population model are described. Section [3] presents the results of the model. In 
Sections |H and E] the results are discussed and summarized. AB magnitudes are used, unless 
otherwise stated, and h = 0.7, Q\ = 1 — Qm = 0.7 is assumed as necessary. 



2. The Model 

2.1. Methodology 

In order to study the average properties of AGN host galaxy stellar populations, a 
model is created which combines knowledge of AGN space density and evolution from deep 
X-ray surveys and CXRB population synthesis models with stellar population studies. A 
model SED is computed which includes the emission from the AGN, the host galaxy stellar 
population, and ongoing dust enshrouded star formation. This AGN and host SED covers 
the wide spectral range of hard X-ray to mid-IR. The calculation of the AGN model SED is 
discussed in Section 12.31 and the details of the host stellar population model SED calculation 
ar e discussed in Sect ion 12.41 To account for dust enshrouded star formation, the templates 



of iRieke et al.l (120091 ) are used, which extend from 5 /mi to 30 cm, and are determined by 
averaging empirical SEDs of local purely star forming galaxies. While dust enshrouded star 
formation may contribute to the galaxy SED at wavelengths shorter than 5 /zm, generally 
the near IR AGN a nd ho st SEDs are dominated by either emission from the AGN or the 



Bruzual fc Charlotl (120031 ) portion of the host SED. Examples of the rest frame AGN and host 



SEDs are shown in Figure [H where the dot-dashed lines sho w AGN SED s , the dashed lines 



show stellar population SEDs, and the dotted lines show the IRieke et al.l ( 120091 ) templates 



A Gaussian stellar mass (M*) distribution was used with M™ 71 = 10 9,5 M Q and M™ ax 
10 12 M . The average M* was set at 10 10,9 M with a standard deviation of 0.4 dex, in 



agreem ent with the sample of X-ray selected, z ~ 1 AGN hosts presented by iPierce et al. 



(120101 ) . It is assumed that the mass-to-light ratio of the host galaxy is constant within 
each waveband and thus the host galaxy flux scales linearly with the host M*. Therefore, 
the stellar population SED is calculated for M* = 1 M Q and then re-normalized to M* = 
10 9 ' 5 - 10 12 M before the stellar population and AGN SEDs are combined. Creating host 
SEDs at several masses in the AGN host mass spectrum allows the consideration of AGN 
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and host colors for hosts of various masses. When considering space densities, however, it 
is important to consider not only the range of AGN host masses but also the probability 
that an AGN host be a certain M*. The space density calculations presented are therefore 
weighted sums across the mass spectrum of AGN hosts where the weight is determined by 
the Gaussian M* distribution described above. Specifically, when calculating luminosity 
functions and number counts the weighted sum over the Gaussian M* distribution is used; 
however, when considering flux-flux plots, only the range of M* is used. In agreement 
with current observational results, the M* distribution used does not evolve with redshift 



flBluck et al.ll201lh 



A variety of broad band data from the optical through mid-IR is used to constrain 
the host models. These model constraints are described in Section 12.51 and include mid-IR 
number counts, the J band AGN and host space density as a function of redshift and the 
absolute J band magnitude, Mj, and the type 1 AGN and host B ban d luminosity function. 
We begin by assuming that the unified model (e.g., lAntonuccil Il993l ) holds and therefore, 
since all AGN are essentially the same system viewed along different lines of sight, there is 
only one class of AGN host galaxies. 



2.2. CXRB Synthesis Model 

It is known that the CXRB encodes the accretion hist ory of SMBHs and as su ch is a 
powerful tool for scientific inquiry into accretion processes (IFabian fc Barconslll992l ). The 
CXRB synthesis model framework can be used in wavelength regions outside of the X-ray 



spectral region by modeling the SEDs of AGIN 
CXRB synthesis modeling (e.g. 



Risaliti et al- 



and their host galaxies. Thus the advances in 



19991 : iBallantvne et al.ll2006al : iGilli et al.ll2007l : 



Treister et al.ll2009bt Draper k, Ballantynd 120091 ) can be used to further the understanding 



of AGN at all wavelengths. Particularly, CXRB synthesis models provide constraints on the 
distribution of N H , the fraction of type 2 AGN, f 2 , and the fraction of CT AGN, f CT . Here 
a simple N H distribution is assumed where type 1 AGN are evenly distributed among log N H 
= 20.0, 20.5, 21.0, 21.5, type 2 AGN are evenly distributed among log A^ = 22.0, 22.5, 23.0 
23.5, an d CT AGN are evenly distributed among log Nh = 24.0, 24.5, 25.0. IBallantvne et al. 



( ]2006aj ) showed that this Nh distribution is consistent with CXRB observation s and provides 



result s which are negligibly diff erent from the observed N H distribution of iRisaliti et al. 



(119991 ) ( iBallantvne et al.l l2006al Jbl). The shap e of the CXRB combin ed with deep X-ray 



surveys has also been used to constrain / 2 
redshift, ft oc (1 -+- z) a (logLx)~ 6 where a 
normalization of ji is set such that at z - 



e.g 



Ball antvne et al.ll2006af). At a given Lx and 
0.4 IBallantvne et al.ll2006ah and b = 4.7. The 



and logLx = 41.5, the type 2 to type 1 ratio 



-6- 



is 4:1 (see Section 2.2 of iDraper &: Ballantyndl2009l ) 



Both the original and the composite model of CT AGN evolution of lDraper fc Ballantyne 



(120 lOl ) are considered. The original, non-evolving model assumes that CT AGN are a sim- 
ple extension of t he Compton thin type 2 AGN population and that, when assuming the 
Ueda et al.l ( 120031 ) hard X-ray luminosity function, ~44% of all obscured AGN are CT. In this 
model there is a population of CT AGN which are nearly as common as Compton thin type 2 
AGN and which are in all ways similar to Compton thin type 2 AGN except for the presence 
of more obscuration along our line of sight. Contrastingly, the composite model assumes that 
CT AGN are a population of AGN distinct from the Compton thin type 2 AGN population. 
Simulations show that gas rich galaxy mergers will cause dust and gas to be funneled into the 
nucle ar region of the galax y, triggering star formation and accretion onto the central SMBH 



(e.g., iHopkins et al.l 120061 ). Due to the l arge reservoir of material, the SMBH will accrete 
very rapidly and be very highly obscured (JFabianlll999l ; Fabian et al.ll2008l 12009). Th us it is 
expected that CT AGN are high Eddington ratio sources. Draper fc Ballantynd (120101 ) found 
that if all CT AGN are rapidly accreting sources, than the local space density of CT AGN 
is under-predicted, and thus a population of low Eddington ratio CT AGN are necessary 
to explain the observed space density of CT AGN. These low accretion rate CT AGN are 
likely obscured by molecular clouds within the host bulge. Observation al evidence of this low 
luminosity CT AGN population h as been found i n the local universe (jTerashima fc Wilson 
20031 ; iGonzalez-Martm et al.ll2009l ). Furthermore, IDraper &: Ballantynd ( 20101 ) showed that 
if CT AGN have moderate Eddington ratios, the space density of CT AGN with Lx > 10 43 
erg s _1 at z > is greatly over-predicted. Therefore, the evolving model of CT AGN has 
an Eddington ratio dependent for with ~86% of all AGN which are accreting at >90% 
of their Eddington rate being CT and ~60% of all AGN which are accreting at <1% of 
their Eddington rate being CT. In this model CT AGN are a distinct population of AGN 
associated with specific stages of AGN evolution and fcr evolves with both L x , the AGN 
2-10 keV luminosity, and redshift. Both the evolving model and the non-evolving model are 
consistent with the CXRB and the local space density of CT AGN with Lx > 10 43 erg s _1 



(IDraper fc Ballantynd 120101 ) . 



2.3. AGN SEDs 



In order to compute the m odel AGN SEDs, the photoionization code Cloudy ver- 
sion C08.00 (IFerland et al.lll9981) is used, following the same procedure as in Section 2 of 



Draper fc Ballantynd (1201 ll ). These SEDs cover the wavelength range of very hard X-ray 
through far-IR and include the transmitted AGN emission, the diffuse emission emitted along 
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the line of sight by the obscuring material around the AGN, and the reflected emission off 
the inner face of the o bscuring cloud. The inner ra dius of the obscuring material is assumed 
to be ~10 pc. As in iDraper fc Ballantynd ( 120111 ). the obscuring clouds are assigned neu- 
tral hydrogen densities, nn, in agreement with typical molecular clouds; the Compton thin 
clouds have Uh = 10 4 cm -3 and the Compton thick clouds have Uh = 10 6 cm -3 . 



2.4. Host Galaxy Stellar Population Model 



The host galaxy stellar populations are modeled using GALAXEV I 1 ] fjBruzual fc Chariot 



20031 ). This model calculates stellar population SEDs over the wavelength range of 91A- 
160 /im. Using an isochrone synthesis technique, GALAXEV can evolve stellar populations 
with three stellar evolution library options. We assume the stellar evolution models of the 
'Padova 1994' library and the stellar spectra of the STELIB/BaSeL 3.1 semi-empirical library. 
GALXEV models can be compute d for six differe n t met allicities in the range Z = O.OO5Z 
- 2.5Z . Both the ISalpeterl (119551 ) and IChabrierl (120031 ) initial ma ss functions (IMFs) are 



available in the GALA XEV co de. Here we ass ume so lar metallicity ( 



Silverman et al.ll2009i ) and the IChabrierl ( 20031 ) IMF ([Silverman et al 



^a uffmann et al.l 120071 ; 



20091 ). The sensitivity 



of the results to these assumptions is assessed in Section [31 GALAXEV also allows for 
a variety of star formation histories, including constant star formation rate, instantaneous 
bursts, and expon entially declining star f ormation rate. An ex ponentially declining star 
formation history (jKauffmann et al.l 120071 ; ISilverman et al.l 120091 ) with an e- folding time of 
0.5 Gyr is used here. This star formation history allows for a simple parametrization of 
the average age of the stellar population while not requiring the entire stellar population 
to be a single age. GALAXEV does not include re-radiation of energy absorbed by dust; 
however, the majority of this re-emission by cool dust occurs at wavelengths longer than 
those considered here. 

The stellar population SEDs are computed in four different redshift bins: z < 1, 1 < z 
< 2, 2 < z < 3, and 3 < z < 5, where z is the redshift. We assume that the host galaxy has 
a dominate stellar population (DSP) which is slightly younger than the average age of the 
universe in each redshift bin, and M* is given by the Gaussian distribution described in Sec 
tion 12 .11 A younger stellar population (YSP) which accounts for ~ 0.1M* (jKauffmann et al. 



2007 



Shi et al.ll2009 ) is also included. The age of t he YSP is selected such that -D n (4000A) 



1.5 (IKauffmann et al.112003 : ISilverman et aUkooa ). using D n (4000A) = fjtc/fuc, where fsc 
is the flux in the 4000-4100 A wavelength range and fsc is the flux in the 3850-3950 A wave- 



available at http://www cida.ve/ bruzual/bc2003 



length range ( Balogh et al.lll999l ). For the stellar population models used here, -D n (4000A) 
« 1.5 corresponds to a stellar population age of ~2 Gyr. The age of the YSP and DSP 
in each redshift bin is summarized in Table [TJ At z > 2, the older and younger stellar 
populations are the same age and thus at 2 < z < 3 the entire stellar population is 2 Gyr 
old and at 3 < z < 5 the entire stellar population is 1 Gyr old. Extinction due to dust, as 
described by E(B — V), is allowed to vary between redshift bins, but not within redshift 
bins. While in some objects the AGN emission is likely to be extin cted by extended dust 
structures within the host galaxy (e.g.. lMartinez-Sansigre et al.ll2010l ). this is not true for all 
objects. As this study focuses on modeling the average properties of AGN and their hosts, 
a detailed modeling of the complex geometry of obscuring material in individual sources is 
beyond the scope of this work. Thus, it is assumed here that all of the extinction suffered 
by the AGN emission occurs within a few tens of parsecs of the central engine. The dust 
which extincts the stellar population is not allowed to further extinct the AGN emission. 



2.5. Model Constraints 



A variety of data in a broad range of spectral regions was used to constrain the stellar 
population models. The number counts are computed at 3.6, 5.7, 8.0 and 24 /im and com- 
pared to the AGN observed by Spitzer Space Observatory in the GOODS field s with / : 
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> 



1 x 10 erg s cm , where /2-8 is the 2-8 keV flux (ITreister et al.l 120061 ). This X-ray 
flux limit is taken into consideration in the number counts calculated here, thus the number 
counts presented are computed for the same population probed by the observational data 
points. These number counts are dominated by type 2 AGN and their hosts and therefore 
are primarily useful to constrain the type 2 AGN host population. The number of sources 
at wavelength A per square degree with flux greater than S, N x (> S), is found by 



iV A (> S) 



K d s ^c 



x 



d$(L x ,z) 



<$ 



max (log L™ n ,log 



L s x) d\ogL x (1 + 3)2^(1 + 3)3 + n A ]V: 



;d log L x dz, 



T 



w here K^ 9 = 3.05 x 10 4 deg 2 sr 1 , dQ/dlogLx is the hard X-ray luminosity function 



of lUeda et al.l (120031 ) or, when r eferring to the evolving mod el, is the evolving Eddington 
ratio space density calculated by iDraper fc Ballantynd ( I2010l ). in units of Mpc -3 , di is the 
luminosity distance, and log L x is the 2-10 keV rest-frame luminosity which corresponds to 
the observed- frame flux S at redshift z. 

As the number counts are dominated by type 2 AGN and their hosts, a separate data 
set must be used to constrain the type 1 AGN host stellar populations. The most obvious 
choice is the optical luminosity function for type 1 Seyferts and QSOs and their hosts. 
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Using the X-ray luminosity function, the type 1 AGN and host B band luminosity function, 
d& B /d(magB) can be calculated using the following equation: 



d$ B d* d(logL x ) 



d(magB) d(\ogLx) d^mags) 

This luminosity function is considered in three redshift bins — z < 0.4, 1.0 < z < 1.55, and 
1.55< z < 2.1. 

Further constraints were considered in order to rule out possible degeneracies between 
stellar population age and extinction due to absorption by dust. The first of these additional 
constraints was the AGN and host J band space density as a function of redshift in several 
Mj bins. First, the AGN and host J band luminosity function, d&j/d(magj), must be 
calculated. This is done in the same manner as d<&B/d(magB) above. Therefore, 

d^j{Mj,z) _ A ^ rf$ d(\ogL x ) 
d(magj) d(logL x ) d(magj) ' 

The normalization constant A depends on the AGN spectral type being considered. For type 
1 AGN, A = (1.0 - / 2 ), while for type 2 AGN, A = f 2 , and for CT A GN, A = fr r. For 



comparison with the J band AGN and host space density measured by lAssef et al.l ( 1201 ll ) , 
d&j/d(magj) is then binned into six Mj bins. 

Additionally, the optical colors were considered in the form of the U-B versus Mg color- 
magnitude diagram (CMD). Also, it is well do cumented that AGN and their hosts tend 



to have X-ray to o ptical flux ratios of Oil (e.g. lAlexander et al.l 12001c iBrandt et al.l 12001 



Rovilos et al.l I201CH ). with CT AGN hosts generally falling below this ratio. Here the ratio 
between the soft X-ray flux, /0.5-2, and the R band flux, /#, log(/o.5-2//i?) = 0±1 is used 
to ensure that the absolute optical fluxes are in agreement with observations. 



2.6. Procedure 

With all the necessary ingredients in place, we begin by assuming that the unified model 
holds, and therefore, on average, type 1, type 2, and CT AGN hosts are identical. The AGN 
space density, Nh distribution, type 1/type 2 AGN ratio, and for are set by the CXRB 
model used. The age of the stellar populations is assigned such that the dominant stellar 
population age is slightly less than the mean Hubble time in each redshift bin and a younger 
stellar population is assigned the age corresponding to the average observed -D n (4000 A) for 
AGN host galaxies. The stellar population ages are summarized in Table [lj The host M* 
distribution is Gaussian within the observed AGN host M* range. The only free parameters 
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are the host galaxy dust enshrouded star formation rate and the extinction due to dust. These 
two free parameters are used to fit the near and mid-IR AGN and host n umber counts. The 
dust enshrouded star formation rate is determined by cycling through the lRieke et al.l (120091 ) 
templates in order of increasing infrared luminosity in steps corresponding to a change in 
SFR of ~1 M Q yr _1 , until the 24 /mi number counts are over-estimated at all flux levels. The 
template of the highest luminosity to not over-estim ate the mid-IR number cou nts is selected 
as the best fit template. The dust extinction in the iBruzual fc Charlotl (120031 ) models is set 
using the total effective V band optical depth obscuring young stars, Ty- SEDs are computed 
with Ty = 0.0, 1.0, 2.0, 5.0, 7.5, 10, 15, 20, and 25. This corresponds to E(B - V) = 0.0 - 
~1.0, with the exact values of E(B — V) depending on the age of the stellar population. The 
E(B — V) with a reduced x 2 , xled^ cl° ses t to 1.0, with respect to the 35 number counts data 
points at 3.6 /mi, 5.7 /mi, 8.0 /mi, and 24 /mi, is selected as the best fit model. The resulting 
model is then compared to the suite of optical and near IR observations described in Section 
12.51 in order to evaluate the appropriateness of the model fit to observations. In this manner 
the average AGN host galaxy properties are elucidated without the complications which 
arise from object-to-object variability. Next, the assumption of the unified model is tested 
by considering the specific constraints for different spectral types. Multiple star formation 
scenarios and the evolution of f% are also considered. 



3. Results 



3.1. Unified Model of AGN Hosts 



In this section, it is assumed that the unified model holds and therefore type 1, type 2, 
and CT AGN have, on average, identical host galaxies. First, the maximum dust enshrouded 
star formation rate (SFR) allowed by the near and mid-IR number counts for X-r a y selec ted 
AGN observed by Spitzer in the GOODS fields, as reported by iTreister et al.l (120061 ) . is 
determined. Figure [1] shows that the mid-IR emission is due primarily to the AGN and dust 
enshrouded star formation; therefore, fitting the 24 /mi number counts can be used to set an 
upper limit on the average AGN host dust enshrouded SFR. In order to not over-estimate 
the AGN and host mid-IR number counts, the avera ge AGN host must have a SFR < 2 
M yr _1 . This SFR is consistent with the findings of iBallantyne et al.l (J2006bf ). who found 
that an AGN host SFR m 1 M Q yr _1 provides a good fit to Spitzer's measurement of the 
AGN contribution to mid-IR portion of the cosmic infrared background, and is similar to 
the SF R of local normal spiral galaxies (e.g., ILee et al.ll2009l ). According to the iKennicutt 
(119981 ) relation, this SFR corresponds to an infrared star formation luminosity Lj^ rj 10 10 
L - 
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Next, the dust content of the average AGN host galaxy is fixed by fitting the near IR 
AGN and host number counts. It is found that if the average AGN host has E(B — V) ~ 
0.5, the near IR number counts are under-predicted. However, if the average AGN host has 
E(B — V) ~ 0.4, the near and mid-IR number counts are over-predict ed. As it is known 



that higher redshif t galaxies tend to be dustier than local galaxies (e.g., iDunne et al.l 12003 



Santini et al.l 120101 ). it is assumed that z < 1 AGN hosts contain less dust than z > 1 AGN 
hosts. If the average AGN host at z < 1 has E(B — V) ps 0.25 and the average z > 1 AGN 
host has E{B — V) ~ 0.5, the predicted number counts are in good agreement with the near 
and mid-IR number counts observations, with xled = 1-2- If the average z < 1 AGN host 
has E(B — V) < 0.25, the mid-IR number counts are over-predicted. In contrast, the model 
constraints are not very sensitive to the age of the older, dominant stellar population, which 
can be changed by ~ 1 Gyr with only minor effects. The age of the younger stellar mass 
population is constrained to within ~0.5 Gyr. 

This unified model of AGN hosts is also in good agreement with the other model con- 
straints. The J band AGN and host space density is in decent agreemen t with obser v ations , 
as shown in Figure [2J The data points in Figure [2] sh ow the resul t s of lAssef et al.l ( 1201 ll ). 
At z < 0.75 the IR and X-ray selected AGN sample of lAssef et al.l (1201 ll ) was chosen using 
a specific selection criterion for optically extended sources. As the models used here do not 
account for the spatial extent of the AGN host galaxies, this exact selection criterion can- 
not be replicated. This gives rise to an obvious discrepancy betw een the mod e l pred ictions 
presented here and the observed J band space density reported by lAssef et al.l (j201lf ) . Thus 
at z < 0.75 we over-predict the AGN and host space density in magnitude bins of medium 
brightness and under-predict the fainter magnitude bins. Type 2 AGN hosts dominate the 
J band space density at Mj > -26 for z > 0.75, with type 1 AGN hosts dominating the 
brightest magnitude bins. 



According to Figure 6 of iPierce et al.l (]2010h . type 1 and type 2 AGN tend to be found 
in the 0.55 < U-B < 1.4 and -23.5 < M B < -19 region of the CMD. This region of the CMD 
is well populated by the model adopted here as many combinations of AGN Lx and host 
galaxy M* fit this criteria. Furthermore, the AGN and hosts fit the expected log(/o.5-2//i?) 
ratio, as shown by the green triangles in Figure El Thus the unified model of AGN hosts 
explains the average trends of AGN and host observations in the optical through mid-IR 
spectral regions. 

In order to understand how well constrained these findings are, the same procedure is 
foll owed using a su persolar metallicity, Z = 2.5Z Q , a subsolar metallicity, Z = 0.2Z & , and 
the ISalpeterl ( 119551 ) IMF. In all of these scenarios, it is found that the stellar ages described 
in Section |2~41 are able to fit the model constraints well; however different E(B — V) values 
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are necessary. The columns labeled " Unified/type 2 AGN ho st" in Table [2] sum marize the 



E(B — V) values for the different models considered. As the iRieke et al.l ( 120091 ) templates 
are independent of the GALAXEV stellar population models, the average AGN host dust 
enshrouded SFR is not affected by changing the metallicity or IMF of the host galaxy stellar 
population. 

Overall, the near and mid-IR AGN and host number counts, as well as the J band 
space density, optical colors, and X-ray to optical ratio of AGN and their hosts are in good 
agreement with observations if type 1, type 2, and CT AGN have, on average, similar host 
galaxies. Also, the near and mid-IR number counts require AGN hosts at z > 1 be dustier 
than AGN hosts at z < 1. We now investigate if this finding holds when spectral type specific 
model constraints are considered. 



3.2. Type 1 AGN Hosts 

In order to study the host galaxies of type 1 AGN, we begin by calculating the type 1 
AGN contribution to the type 1 AGN and host B band luminosity function. The predicted 
AGN only B band luminosity function is shown as the dot-dashed lines in Figure HI where 
black lines and data refer to z < 0.4, blue lines and data refer to 1.0 < z < 1.55, and red 
lines and data refer to 1.55 < z < 2.1. The Xred takes into account the 50 B band type 1 
AGN and host B band luminosity function data points shown in Figure HI When only the 
AGN contribution to the B band luminosity function is considered, xled = 14- It is clear that 
except for at the very brightest magnitudes, a contribution from the host galaxy is necessary 
in order to fit the observed type 1 AGN and host B band luminosity function. 

If the same stellar population described in Section 13.11 is used to calculate the type 1 
AGN and host B band luminosity function, the predicted z > 1 luminosity functions are 
not in agreement with observations. The unified model gives a Xred = 13 fit to the B band 
luminosity function. This is shown by the dashed lines in Figure SJ For the z < 0.4 type 1 
AGN and hosts, the lower limit E(B — V) of the unified AGN host model B band luminosity 
function (solid black line in Figure H]) is in decent agreement with observations. As the type 
2 AGN hosts are constrained by the IR AGN and host number counts as described in Section 
13. 1[ this suggests that at z < 1, AGN hosts can be described by the unified AGN host model, 
but at z > 1 the unified AGN host model cannot properly predict the type 1 AGN and host 
B band luminosity function. Thus at z > 1, type 1 and type 2 AGN must be hosted by two 
distinct galaxy populations. The lowest Xred "with respect to the B band luminosity function 
is found when type 1 AGN hosts have E(B — V) ~ 0.05 at z > 1. For this model xled = 4-3, 
which is not formally a good fit to the data points, but is dominated by the over-prediction 
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of the two faintest ICroom et al.l (120041 ) data points at both 1.0 < z < 1.55 and 1.55 < z < 
2.1. If these points are not included Xred = 2-5. As these data points refer to observations 
of faint AGN at high redshift, it is likely that the observational sample is not complete at 
these fluxes. Thus, at z < 1 type 1 and type 2 AGN hosts are similarly dusty. However, at 
z > 1, type 1 AGN hosts are less dusty than type 2 AGN hosts. Also, on average, it appears 
that type 1 and type 2 AGN hosts have similarly aged stellar populations, at least at z < 2, 
where the model can be well constrained. 

The type 1 AGN host stellar population described above is in good agreement with the 
other model constraints considered. Figure |5] shows the number counts for type 1 AGN and 
their hosts as dotted lines. This model has a Xred = 1-3 with respect to the observed number 
counts. For the unified model, a xled = 1-2 was found with respect to the IR number counts; 
thus reducing the dust extinction in type 1 AGN hosts at z > 1 has little effect on the total 
IR number counts. Type 1 AGN and hosts make a considerable contribution to the shorter 
wavelength number counts, but at mid-IR wavelengths they make a minimal contribution. 
Again shown as dotted lines, Figure [6] demonstrates that type 1 AGN and hosts dominate 
the J band space density for Mj < -26 and make a significant contribution in the -26 < Mj 
< -24 magnitude bin. However, at fainter absolute magnitudes, the type 1 AGN hosts make 
only a nominal contribution to the J band AGN and host space density. 

The region of the CMD where type 1 AGN and the i r host s tend to be found, as shown by 



the upside down triangles in Figure 6a of IPierce et al.l ( 120101 ) . is -23.5 < Mb < -19 and 0.55 



< U-B < 1.4. This region of the CMD is well populated by a variety of AGN Lx and host 
M* model combinations. The majority of type 1 AGN and host model SEDs used here are in 
agreement with the log(/o.5-2//R) ratio, as exhibited by the blue circles in Figure[3j At fainter 
soft X-ray fluxes it appears that the type 1 AGN and hosts fall below the expected ratio. 
These sources are located at z > 2.0 where the observed R band flux measures rest frame UV 
emission. Due to the observed frame R band filter redshifting out of optical wavelengths, it 
is unclear whether the log(/o.5-2//ij) = ± 1 relation should hold at high redshift. These 
high redshift type 1 AGN and hosts can be brought into agreement with log(/o.5-2//i?) — 
± 1 by increasing the extinction due to dust in the host stellar population or by increasing 
the age of the host stellar population. Thus if the locally observed log(/o.5-2//ij) ratio holds 
at high redshift, this suggests that at the peak of quasar activity, type 1 AGN hosts either 
had older stellar populations than type 2 AGN hosts at the same redshift or were dustier 
than 1.0 < z < 2.0 type 1 AGN hosts. If the stellar populations of type 1 AGN hosts at 
z > 2 are older than the stellar populations of type 2 AGN hosts at the same epoch, this 
implies an evolutionary scenario where recently triggered AGN are obscured and then blow 
out the obscuring gas and dust in order to reveal an unobscured AGN and an older host. It 
is also possible that, since the fraction of gas rich galaxies at z > 2 is larger than at 1 < z 
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< 2 fjDahlen et al.ll2007l ). type 1 AGN hosts located at z > 2 contain more dust than type 1 
AGN hosts located at 1 < z < 2. Both of these possibilities are consistent with quasars being 
triggered by major mergers, which is expected to be the dominate quasar fueling mechanism 
at 2 < z < 3. 



When considering the subsolar metallicity, supersolar metallicity and ISalpeterl (119551 ) 
IMF models, the values for E(B — V) which supply the best fit for the type 1 AGN and 
host B band luminosity function change, but the general trend that type 1 and type 2 AGN 
host stellar populations are similar at z < 1 and different at z > 1 remains. The E(B — V) 
values for type 1 AGN hosts in the different models are shown in the columns labeled "Type 
1 AGN host" in Table [2j In all models the type 1 AGN hosts are similar in dust content to 
the type 2 AGN at z < 1 and the type 1 AGN hosts are less dusty than the type 2 AGN 
hosts at z > 1. Thus, regardless of the metallicity or IMF, type 1 and type 2 AGN hosts are 
similar at z < 1 and intrinsically different at z > 1. 



3.3. Star Formation in AGN Hosts 



The emission due to obscured star formation is taken into account using the iRieke et al. 



( 120091 ) star formation templates. The mid-IR AGN and host number counts suggest that 
the average AGN ho st has a dust ob scured SFR « 2 M Q yr" 1 , the same SFR as normal local 
spiral galaxies (e.g., [Lee et al.ll2009l ). Ho wever, recent studies suggest that a large r, but still 
modest, AGN host SFR is expected (e.g. jLutz et al.ll2010l ; Melbourne et al.ll2011). Further 



more , studies also find t hat the AGN host SFR tends to increase with redshift ( 



2010). AGN lumino sity ([Thompson et all 120091 ). or both (jSerjeant fc Hatziminaoglou 



Lutz 



et al. 



2009 



Serjeant et a 



of 



20101) . As several AGN fueling mechanisms, such as the starburst disk model 



Ballantynd ( 20081 ). require processes related to star formation, it is important to consider if 
it is possible for AGN hosts to have average SFR > 2 M yr _1 . Thus several SFR evolutions 



are considered here. 



The lWilman et al.l (120101 ) SFR evolution depends on both AGN Lx and redshift, finding 

SFRoc v / Wl0 43 (1.0 + 2) L6 . (4) 



When the lWilman et al.l (120101 ) SFR evolution is used, the average z < 1 SFR must still be < 
2 M yr _1 in order to not over-predict the mid-IR number counts. However samples of AGN 
hosts with aver age SFRs an order of magnitude higher than this (~ 18-41 M ^ yr _1 ) have 



been observed ( ILutz et al.l 12010c Melbourne et al.l 12011c ISeymour et al.l 120111 ). In order to 
explain these observations, a population of enhanced star formation sources are considered. 
The normalization factor for the proportionality in equation [4] is set such that the average 
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SFR for AGN hosts at z < 1 is ~20 M yr _1 . In order to not over-estimate the faint end of 
the 24 fim number counts, enhanced star formation sources can account for at most ~15% 
of the AGN population. 

Si milar fractions of enh anced star formation sources are found using other SFR evolu- 



tions. ISerjeant et al.l ( 120101 ) used Herschel Space Observatory observations of type 1 SDSS 
selected quasars to study the evolution of SFR in quasar hosts, finding that AGN host SFR 
displays strong luminosity dependent evolution with redshift. This leads to a SFR evolution 
of the form 

SFR oc (1.0 + z) Q , (5) 

where a « -1.91 ab - 42, where Iab is the absolute I band magnitude of the quasar. If the 
average z < 1 SFR is set at ~20 M yr _1 , enhanced star formation sources can at most be 
~5% of the AGN population. 

If AGN host SFR does not evolve with luminosity, and instead only evolves with redshift, 
then the average AGN host SFR can be considera bly higher. If the average AGN host SFR 
redshift evolution found by ISerjeant et al.l (120 10[ ) for moderate luminosity AGN, SFR ex 
(1.0 + z) 2 - 3 , is applied to all AGN hosts, an average z < 1 SFR « 16 M yr _1 for both type 
1 and type 2 AGN is in good agreement with the mid-IR AGN and host number counts, 
with xled = 1-3- This avera ge SFR co r respon ds to SFR(,z=0.0) = 0.5 M Q yr _1 , which is 



in excellent agreement with iKim et al.l (120061 ) who, using a sample of local type 1 SDSS 
quasars, found that the average local AGN host SFR fh 0.5 M Q yr _1 . Furthermore, in the 
redshift only SFR evolution, for an average z < 1 SFR s=s 20 M yr _1 , ~80% of AGN hosts 
can be enhanced star formation sources. 

Thus, if the AGN host SFR evolution is dependent on luminosity, either through lumi- 
nosity evolution or luminosity dependent redshift evolution, enhanced star formation sources 
are ~5-15% of the AGN population, but if AGN host SFR only evolves with redshift, the 
average AGN host SFR can be up to ~16 M yr _1 . This finding is fully consistent with 
the 15 /im AGN luminosity function at z ~ 0.7 observed by iFu et al.l (120101 ) . Furthermore, 
Lutz et al.l ( 120101 ) find that for moderate luminosity AGN, the SFR does not evolve with Lx 
and that only at the highest quasar luminosities does the AGN host SFR seem to depend on 
the AGN luminosity. Therefore, the fraction of enhanced star formation sources found here 
is a lower limit to the true fraction of enhanced star formation sources. 
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3.4. CT AGN Hosts 

Observational and theoretical evidence suggests that CT levels of obscuration of an 



AGN may be due to an evo 
phase of rapid growth (e .g 



2008; Fabian et al. 2009 



utionary stage whe r e the SMBH and host bulge a r e both in a 



Sanders et al 



1988; 



Fabian 



1999 



Page et al.l 120041 ; iBallantvne 



Draper &: Ballantynd |2010| ; iNardini fc Risalitil l201ll ) . According 



to galaxy merger simulation s, gas rich mergers will ignite a burst of star formation and 
rapid black hole growth (e.g. iHopkins et al.l 120061 ) . In this scenario, it is expected that CT 
AGN hosts would be a subset of the ultraluminous infrared galaxy (ULIRG) population, 
characterized by Lj^ > 10 12 L Q , where Ljr is the 8-1000 /zm luminosity. In order to test 
this scenario, we attempt to model the CT AGN hosts with ste ll ar po pulations similar to 



those found in ULIRGs hosting AGNs. iRodriguez Zaurin et al.l (120101 ) find that ULIRGs 



hosting AGNs have an average stellar population age of ~0.3 Gyr which dominates the 
stellar mass of the galaxy with an average M* « 10 ia8 M and standard deviation of ~0.35 
dex and an average E(B — V) ~ 0.6. This scenario is tested using both the non-evolving 
and evolving models described in Section 12.21 

Due to the observational challenges of identifying CT AGN, there is no data set specific 
to the optical or near infrared properties of CT AGN hosts. Thus the near and mid-IR AGN 
number counts measured from an X -ray flux limited sample of AGN, with / 2 -8 ^ 1 x 10~ 16 
erg s _1 cm -2 (ITreister et al.l 120061 ). are used to place limits on the average CT AGN host 
galaxy. The X-ray flux limit of the Spitzer GOODS AGN sample is taken into account in 
the number counts calculation . Thus, despite the fact that a large fraction of CT AGN are 
missed by deep X-ray surveys (lHasingerll2008l ). the d ot-dashed lines in Fig ures [5] and [7] show 
the CT AGN which would have been selected by the ITreister et al.l (120061 ) selection criteria. 
The J band space density, X-ray to optical flux ratio, and optical colors are also used to 
ensure the model CT AGN hosts are consistent with observations of the AGN population as 
a whole. 



3.4-1- Non- evolving Model 



The non-evolving model of iDraper fc Ballantynd (120101 ) assumes that CT AGN evolve 
like less obscured type 2 AGN. In order to match the peak of the CXRB at ~30 keV, it is 
required that ~44% of obscured AGN are C T. When the CT AGN hos t stella r populations are 
modeled in agreement with the findings of IRodriguez Zaurin et al.l (120101 ). the IR number 
counts are greatly over-predicted. This suggests that on average, CT AGN in the non- 
evolving model cannot be hosted by ULIRGs. In fact, when using the non-evolving model, 
the average CT AGN host cannot even be a luminous infrared galaxy (LIRG), which is 
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characterized by Lin > 10 11 L . In order for the near and mid-IR number counts to not be 
over-predicted at the faint end, the average non-evolving model CT AGN host must have 
Lj r 5, 10 10 - 75 Lm, co rresponding to a star formation rate of < 10 M Q yr -1 , according to 



the iKennicuttl (119981 ) relation. The youngest the non-evolving model CT AGN host average 
stellar population can be, and not over-predict the mid-IR number counts, is ~ 1 Gyr old 
which requires E(B — V) ~ 1.0. The contribution of the non-evolving model CT AGN and 
their hosts to the near and mid-IR number counts is shown as the dot-dashed line in Figure 

El 

In the non-evolving model, CT AGN dominate the low redshift, faint Mj region of 
the J band space density. Figure [6] shows the CT AGN and host contribution to the J 
band space density as the dot-dashed lines. As expected by the higher dust content, CT 
AGN hosts occupy a region of the CMD which is on average slightly redder and fainter 
than the type 1 and type 2 AGN hosts. The CT AGN tend to lie below the log(/o.5-2//ii) 
ratio, as the majority of the soft X-ray flux of CT sources is absorbed by the CT obscuring 
material. Figure [3] shows the log(/o.5-2//i?) ratio of CT AGN and their hosts as red squares. 
Changing the metallicity or IMF used for the host galaxy stellar population does not change 
these results. 

When the enhanced star formation sources are considered, the CT AGN hosts must be 
similar to the type 2 AGN hosts, regardless of host galaxy stellar population metallicity and 
IMF. Similarly, when the redshift evolution of AGN host SFR is included, the CT AGN 
hosts must be similar to the type 2 AGN hosts. Thus, if for does not evolve with AGN 
Eddington ratio, then CT AGN and their hosts are expected to be a simple extension of the 
less obscured type 2 AGN population. This would require that AGN triggered by mergers 
be a small minority of the quasar population. 



3.4.2. Evolving Model 

The evolving model allows CT AGN to evolve independently of the less obscured type 2 
AGN. Instead, fcr is assumed to be Eddington ratio dependent. In order to fit the peak of 
the CXRB at ~30 keV, the local CT AGN space density, and the z > 1 IR CT AGN space 
density, it is found that ~86% of AGN with Eddington ratios >0.9 are CT, ~60% of AGN 
with Eddington ratios < 0.01 are CT, and ~0% of AGN with intermediate Eddington ratios 
are CT. 

As AGN with weaker ac c retion rates are observed to have older stellar populations 



( IKauffmann et al.l 120031 12004 . 120071 ). the low Eddington ratio AGN are assumed to have 
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the same mass distribution as the type 1 and type 2 AGN but an older stellar population 
and very little star formation. The age of the low Eddington ratio CT AGN host stellar 
populations are summarized in Table [TJ If the low Eddington ratio CT AGN hosts have 
the same E(B — V) as the type 1 and type 2 AGN hosts, the near and mid-IR number 
counts are over-predicted. Increasing the dust in the low Eddington ratio CT AGN hosts 
to E(B — V) ~ 1.0 brings the model predictions into agreement with the faint end of the 
IR number counts, as shown by the blue dot-dashed lines in Figure [7J The larger amount 
of dust in the CT AGN hosts as compared to the type 2 AGN hosts makes sense as only 
galaxies containing a large amount of dust will be able to host an AGN with CT levels of 
obscuration. 

In the evolving model, the high Eddington ratio CT AGN hosts should be galaxies in 
a phase of rapid star formation and black hole growth. Therefore, the high Eddington ratio 
CT AGN hosts should be LIRGs or ULIRGs. Indeed we find that if the stellar populations of 
the high Eddington ratio CT AGN hosts have Lm = 10 12 L from star formation, E(B — V) 
~ 0.6, with a stellar population of age ~0.3 Gyr which domin ates M*, in agreement with th e 



average stellar population of ULIRGs which host an AGN (IRodriguez Zaurin et al.ll2010l ). 
the near and mid-IR number counts predictions are in agreement with deep observations. 
The high Eddington ratio CT AGN and host number counts are shown as the red dot- 
dashed lines in Figure [7J The over-prediction at the bright end of the number counts is due 
to th e incompleteness of the survey at bright fluxes. As the GOODS fields cover only 0.1 



deg 2 (ITreister et al.ll2006l ). GOODS misses bright, rare objects. If the stellar population is 
assumed to be younger than ~0.3 Gyr, the number counts are over-predicted. Similarly if the 
E(B — V) < 0.6 or if the average M* > 10 11 M , the number counts are over-predicted. Thus, 
the stellar population described here is the upper limit for how bright the high Eddington 
ratio CT AGN hosts can be without over-predicting the IR number counts. 

The CT AGN host galaxies are also in good agreement with the other model constraints 
considered here. The CT AGN host contribution to the J band space density is shown as 
the dot-dashed lines in Figure [8], where the red lines refer to the high Eddington ratio CT 
AGN and their hosts and the blue lines refer to the low Eddington ratio CT AGN and 
their hosts. As expected the high Eddington ratio sources dominate the CT AGN and host 
contribution in the brighter magnitude bins while the low Eddington ratio sources dominate 
at the fainter magnitude bins. The optical colors are also in agreement with observations. 
The low Eddington ratio CT AGN hosts are on average redder and fainter than the type 
1 and type 2 AGN hosts, while the high Eddington ratio CT AGN hosts are on average as 
bright or slightly brighter than the type 1 and type 2 AGN hosts but on average a little 
redder than the type 1 and type 2 AGN hosts. As with the non-evolving model, most of the 
CT AGN hosts lie below the empirical average log(/o.5-2//i?) ratio, as expected by the high 
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levels of soft X-ray absorption fundamental to CT AGN. In Figure [9] the /r versus /0.5-2 for 
CT AGN and their hosts are shown as the red squares with the low Eddington ratio objects 
shown as filled red squares and the high Eddington objects shown as open red squares. The 
same result is found when the host galaxy IMF or metallicity is varied. 

As the high Eddington ratio CT AGN are a small fraction of the overall AGN population, 
the evolving model CT AGN hosts are not affected by the consideration of the enhanced 
star formation sources. Thus the CT AGN hosts of the evolving model are in agreement 
with the AGN evolution scenario where major mergers trigger nuclear starbursts and highly 
obscured AGN activity. 



3.5. Evolution of f 2 

The evolution of f 2 , the type 2 AGN fraction, is important for understanding the A GN 
life cycle and how AGN and their host galaxies interact (e.g., iBallantyne et al.l l2006ai ). In 
the unified model, f 2 is the covering factor of the dusty to rus, and thus the evolu tion of 
f 2 shows a fundamental evolution of the torus parameters (IBallantyne et al.l l2006bh. Sev- 



eral studies suggest that f 2 evolves with Lx and possibly also with redshift (e. 



2003; ISimpsonll2005l : IBallantyne et al.ll2006ai : lHasingerll2008l : Fabian et al.ll2009 



g-- 



Ueda et al. 



Winter et al. 



200C 



However, other studies find that it is not necessary for f 2 to evol v e with L x (e.g., 



Lawrence fc Elvisl |2010| ) nor with redshift (e.g., iGandhi fc Fabian! 120031 ; iTreister fc Urry 
20051 ) in order to explain observations. Indeed, we find that the number counts and J band 
space density can be fit with a constant f 2 . For z < 0.4, the type 1 AGN and host B band 
luminosity function can also be satisfactorily fit with a constant f 2 . However, fitting the 
type 1 AGN and host B band luminosity function at z > 1 requires that f 2 is not constant. 
If f 2 is constant, the best fit Xred > 7.0 and the type 1 AGN and host B band luminosity 
function predicts a considerably larger population of B band magnitude, uib ~ 23, 1 < z < 
2 type 1 AGN than observed. Thus, it is found that in order to fit the type 1 AGN and host 
B band luminosity function at z > 1, f 2 must evolve with Lx- 

The argument for evolution of f 2 with redshift is less conclusive. Figure [10] shows the 
type 1 AGN and host B band luminosity function assuming that f 2 does not evolve with 



6.5. When f 2 does evolve with redshift, the best fit x 



red 



redshift, which gives Xred 
4.3. While the Xred f° r both the evolving and the non-evolving f 2 models do not represent 
formal good fits to the observed data points, the model in which f 2 does evolve with redshift 
provides a better fit to the data. In the scenario where f 2 does not evolve with redshift, it is 
necessary that at higher redshift the type 1 AGN host galaxies be dustier than locally. This 
is in contrast to the scenario where f 2 does evolve with redshift and type 1 AGN hosts are 
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less dusty at higher redshift compared to the local population. In order to not over-predict 
the IR number counts, the type 2 AGN hosts must also be dustier at higher redshift than 
locally Assuming that / 2 does not evolve with redshift has minimal affect on the J band 
space density, optical colors, and log(/o.5_ 2 //jj) ratio. 

A comparison of Figure 0] and Figure [TOj suggests that / 2 not only evolves with Lx, 
but also with redshift. Also, the model where / 2 evolves with both Lx and redshift provides 
a better fit to the type 1 AGN and host B band luminosity function has measured by 
xled- As neither the evolving nor non-evolving / 2 models provide a formally good fit to the 
observed luminosity function, the possibility that / 2 does not evolve with redshift cannot be 
conclusively ruled out; however, the evolution of the B band type 1 AGN and host luminosity 
function suggests that / 2 does evolve with redshift. In order to fit the type 1 AGN and host 
B band luminosity function f 2 must evolve with Lx- 



3.6. Summary of Results 

Using a variety of optical and near and mid-IR data, we have constrained the stellar 
populations of AGN host galaxies. Table |3] summarizes the average AGN host galaxy for 
different AGN spectral types. It is found that at z < 1 type 1 and type 2 AGN hosts are 
similar, but at z > 1 type 1 AGN hosts are less dusty than type 2 AGN hosts. The majority 
of AGN have an average SFR < 2 M Q yr _1 , however there is evidence of a population of 
enhanced star formation sources which account for > 5-15% of the AGN population and 
has SFR » 20 M yr" 1 . Also, it is found that if CT AGN evolve like type 2 AGN, then CT 
AGN hosts are similar to type 2 AGN hosts and if CT levels of obscuration are indicative of 
specific evolutionary stages in the AGN life cycle, then CT AGN hosts are also in specific 
evolutionary stages. Furthermore, it is shown that / 2 evolves with Lx and / 2 is likely to 
evolve with redshift. 



4. Discussion 

4.1. Hosts of CT AGN 

For both the evolving and non-evolving models, the CT AGN host stellar populations 
suffer from at least as much dust extinction as the type 2 AGN host stellar populations. Ob- 
servations suggest that at least some of the observed extinction of heavily obscur ed AGN may 



be due to extended dust s t ructures or molecular clouds w ithin the host galaxy (IBrand et al. 



20071 ; iPolletta et al.l 120081 ; iMartinez-Sansigre et al.l 120101 ) , so it is expected that the stellar 
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populations of CT AGN will be enshrouded in dust. It is within the limits imposed by the 
model constraints for the high Eddington ratio CT AGN hosts to have an average E(B — V) 
which is in agreement with ULIRGs hosting AGN (JRodriguez Zaurin et al.ll2010l ). 



The AGN evolution scheme, in which mergers trigger large nuclear starbursts and AGN 
activity, claims that the AGN activity will initially be very highly obscured while the black 
hole grows very rapidly. As the black hole grows, the radiation pressure on surrounding dusty 
gas will increase until the AGN feed back blows out the obscuring material and halts the star 



formation in the host nucl ear region (ISanders et al.lll988l ; iPage et al.ll2004J ; iRigopoulou et al. 



20091 ; iHopkins et al.ll2006l ). In this scheme CT AGN should have young stellar populations 
and possibly high levels of on going star formation. The non-evolving model places a lower 
limit on the stellar population age of 1 Gyr. However, when the enhanced star formation 
sources are included, the non-evolving model CT AGN hosts have similar stellar ages to the 
type 2 AGN hosts, in contrast to the expectations of the AGN evolution scheme. For the high 
Eddington ratio CT AGN in the evolving model, the lower limit o n the stellar population 



ageis 0.3 Gyr, the average stellar age of ULIRGs hosting an AGN (JRodriguez Zaurin et al. 



2010|). The low Eddington ratio CT AGN hosts of the evolving model can be of similar age as 
the type 1 and type 2 AGN, but the near and mid-IR predicted number counts are in better 
agreement with the observations if these AGN hosts have slightly older stellar populations 
than the average type 1 and type 2 AGN hosts. It has been demonstrated that galaxies 
hosting AGN with lower [OIII] luminosities have larg er values for D n (4Q00 A) compare d 
to galaxies hosting AGN with high [OIII] luminosities (JKauffmann et all 12003 . 120041 . 120071 ). 
Thus it is expected that AGN with lower Eddington ratios are in hosts with older stellar 
populations, in agreement with the findings of this study. 

The stellar populations of CT AGN hosts for the evolving model and the non-evolving 
model are in agreement with the expectations from the AGN evolution scheme. The non- 
evolving model finds that the average CT AGN host has recently (>1 Gyr ago) undergone 
a large burst of star formation, but that current star formation rates are more modest. The 
evolving model finds that high Eddington ratio CT AGN hosts have recently (>0.3 Gyr ago) 
undergone a large burst of star formation and that current star formation rates may also 
be elevated. If the enhanced star formation sources are included, the non-evolving model 
AGN hosts have stellar populations of similar age as the type 2 AGN, and thus are in better 
agreement with the orientation based unified model than with the AGN evolution scheme. 
The evolving model CT AGN hosts are consistent with the AGN evolution scheme regardless 
of the inclusion of the enhanced star formation sources. 

In summary, the non-evolving model CT AGN hosts are a simple extension of the type 
2 AGN host population while the evolving model CT AGN are consistent with the paradigm 
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where major mergers cause both intense starbursts and AGN activity. In the evolving 
model, low Eddington ratio CT AGN hosts will appear as galaxies with old and dusty stellar 
populations while the high Eddington ratio CT AGN hosts will appear as IR bright starburst 
galaxies. 



4.2. Enhanced Star Formation in AGN Hosts 



It is found that ~5-15% of AGN hosts can be enhanced star formation sources with an 
average SFR « 20 M yr -1 , a factor of 10 higher than the majority of AGN hosts. This 
is expected to be a lower limit of the fraction of enhanced star formation sources as some 
sources may have such highly embedded star formation that the majority of the reprocessed 
emission due to star formation is at wavelengths longer than 24 /im. Observations at longer 
wavelengths, such as in the far-IR with the Herschel or the millimeter/sub-mm regime with 
the Atacama Large Millimeter/submillimeter Array (ALMA), are necessary fo r uncovering 
the evolution of star formation rates in AGN hosts (jDraper fc Ballantyndl2011l ). Even deep 
radio observati ons are a useful t ool in determining the highly embedded star formation rates 
of AG N hosts tea 
AGN, 



Lutz ct al. 



Iantyndl2009f) . Indeed, by stacking sub-mm observations of X-ray selected 



(120 101 ) found that the average AGN host SFR « 30 M Q yr 1 . Furthermore, 



Lutz et al.l (120101 ) found that the AGN host SFR evolves strongly with redshift but with 
evidence of luminosity dependent evolution only for the highest luminosity AGN. In the 
scenario where AGN host SFRs are not luminosity dependent, it is found that more than 
half of AGN hosts can be enhanced star formation sources. 

The fact that the mid-IR AGN and host number counts are over-predicted by an average 
SFR > 2 M Q yr -1 , despite observational evidence that there is a population of AGN with 
an average SFR an order of magnitude higher than this upper limit, suggests that there 
are two populations of AGN. The majority of AGN hosts, at least at z < 1, have SFRs 
similar to local spiral galaxies and >5-15% of AGN hosts have markedly higher SFRs. The 
existence of these two populations of AGN hosts can be interpreted in two complimentary 
ways. The first interpretation is that the AGN hosted by galaxies with enhanced SFRs 
are being fueled by different mechanisms than the lower star formation objects. At z < 
1, the lower star formation objects are likely dying quasars which were triggered by major 
mergers while the enhanced star formation objects are likely Seyfer ts which ar e both fueled 
and obscured by circumnuclear starburst disks, as investigated by iBallantynd (120081 ). The 
other interpretation is that all AGN are fueled by processes related to nuclear starbursts, 
whether those starbursts are triggered by mergers or through secular processes, and on 
average, the nuclear starburst phase overlaps with the active AGN phase for >5-15% of the 
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AGN lifetime. Further exploration of these two populations of AGN hosts is necessary to 
determine the processes which trigger and fuel AGN activity. 



4.3. Methods for Finding CT AGN 

It is well documented that the integrated emission of AGN observed in deep X-ray sur- 
veys is insufficient to account for the intensity of the CXRB at ~30 keV and that the shape 
of the CXRB necessitates that the missing population of AGN be hig hly obscured (e.g., 
Ballantyne et al.ll2006al ; iDraper fc Ballantyne! l2009[ iTreister et al.ll2009bl ). Given the uncer- 
tainties of the normalization of the CXRB and the AGN hard X-ray lumi nosity functions, 
different models predict vastly different numbers of missing CT AGN (e.g.. iGilli et al.ll2007l ; 
Draper fc Ballantyn e 2009; Tre ister et al.ll2009bt) . As these highly obscured AGN are missed 



in deep X-ray surveys ( Hasingerl 120081 ). it is important to consider other methods to identify 



the elusive CT AGN population. 

One possib ility is that the majority of X-ray bright optically inactive galaxies (XBONGs) 
host CT AGN (JFiore et al.ll2008l 120091 : ITreister et al~1l2009al : iRovilos et al.lhoioL XBONGs 
are X-ray sources found in deep surveys which have no optical cou nterparts w i th R < 25.5 
and make up a substantial portion of deep survey X-ray sources (jAird et al.ll2010l ). It is 
thought that these sources are either heavily obscured AGN an d/or high redshift quasars 
( lAlexander et al.l l2001t iMainieri et al.l l2005t IRovilos et al.ll2010l ). Considering the sources 
above the dotted horizontal line in Figures [3] and [9] suggests that a small fraction of optically 
faint X-ray sources may be CT AGN, but the majority of this population is z > 1 type 2 AGN. 
Thus it is unlikely that XBONGs host the majority of the missing CT AGN population. 

Another method used to identify CT AGN candidates is to search for infrared bright 
sources whic h have an infrared excess, usually defined by /24/f/? ->, 1000, where fi\ is the 24 



jUmflux fe.g.lPolletta et a 



20061 : lAlexander et al.ll2008l : iDonlev et al-lbOOa iFiore et al.ll2008 



20091 ; ITreister et al.ll2009al ) . Some concern has been raised that this method will also pick ou t 
lower redshift type 2 AGN masquerading as high redshift CT AGN ( IGeorgakakis et al.ll2010l ). 
We therefore investigate the population of AGN selected by the infrared excess criteria. 



Fiore et al.l ( 120081 ) suggest that using the criteria /24//R > 1000 and R - K > 4.5 



selects a distinct class of sources, the majority of which are CT AGN. However, according 
to both the evolv ing and non-evolving model CT AGN and host SEDs developed here, the 
Fiore et al.l (120081 ) criteria selects low to moderately X-ray bright AGN at moderate to high 
redshift. The vast majority o f the selected AGN are obscured, but a significant fr action are 
still Compton thin. Similarly, IDonlev et al.l ( 120081 ) find that the lFiore et al.l ( 120081 ) criteria is 
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l ikely to select dusty star forming templates and low X-ray flux AGN. However, iDonley et al. 



( 120081 ) find no evidence that the selected AGN are CT. 

A m ore effective method of identifying CT AGN candidates is based on the selection 



criteria of iPolletta et al.l (120081 ): f^j fa > 1000 and / 2 4 > 1.0 mJy. Furthermore, it appears 



that samples with a lower / 24 limit, such as f 2 4 > 700 /xJy or even / 24 > 550 /iJy, also contain 



a large fraction of highly obscured AGN ( iDonley et al.ll2010L and references therein). Fainter 



samples of infrared bright g alaxies are found to be predominately powered by star formation 



rather than by AGN (e.g., iPope et al.ll2008l ). Figure ITTa shows that representatives of all 
spectral types of AGN can be found with /24//.R > 1000 and / 2 4 > 550 /iJy. These sources 
will be located at all redshifts. However, the majority of type 1 and type 2 AGN are bright in 
the soft X-ray. For both the evolving and non-evolving models, the vast majority of infrared 
excess AGN with / 24 > 550 /iJy and fo.5- 2 < 10~ 15 erg cm -2 s _1 are CT AGN, as shown in 
Figure [TTb . Figure [12b, which shows the redshift distribution of the CT AGN with fo.5_ 2 < 
10~ 15 erg cm~ 2 s _1 , illustrates that these X-ray faint CT AGN are located at all redshifts. 
For the evolving model the vast majority of these CT AGN are high Eddington ratio sources. 
This criteria will select a small number of type 2 AGN with enhanced star formation, but 
the majority of sources selected in this manner are indeed CT. 

If an infrared excess AGN has fo.5-2 ^ 10~ 16 erg cm -2 s _1 and / 2 4 > 550 /j,Jy, than the 
AGN is CT and located at z < 2. This is true for both the evolving and non-evolving models. 
An interesting consequence of this is that a size-able population of z > 2 CT AGN should 
have 10~ 15 > /0.5-2 > 10 -16 erg s _1 cm -2 . According to the evolving model, nearly all of 
the /0.5-2 ;$ 10~ 16 erg cm -2 s _1 sources are high Eddington ratio sources and according to 
the non-evolving model these low X-ray flux sources are nearly all enhanced star formation 
sources. It is therefore expected that X-ray stacking of IR bright sources with /24//J? ^ 1000 
and low soft X-ray flux will yield a large fraction of the z < 2 high luminosity CT AGN 
population. This shows that combining observations in multiple spectral regions, such as 
mid-IR, optical, and X-ray observations, is the most efficient way of identifying CT AGN 
candidates. 

Identifying and characterizing the elusive CT AGN population is a necessary part of 
understanding the history of accretion and galaxy evolution. However, due to the obser- 
vational challenges of studying highly obscured sources, multi-wavelength investigations are 
necessary to identify and understand the nature of CT AGN. 
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4.4. Implications for the Unified Model and AGN Fueling Mechanisms 

Several recent studies suggest that there are t wo distinct processes which lead to AGN 



activity, secular evolution and merger event s (e.g., iBallantyne et al.ll2006al ; iHasingerl 12008 

Hopkins fc Hernquistl 120091 ; iLutz et al.ll2010l ). In the latter paradigm, quasar activity is ac- 

tivat ed by galaxy mergers which cause gas and dust to be tunneled into the nuclear region 



Sanders et al.lll988t lFabianlll999t iPage et al.ll2004j ; iHopkins et al.ll2006l ). whereas mod 



(e.g. 

erate luminosity AGN are fueled by gravitational instabilit i es internal to the host galaxy 



or through minor interactions (e.g.. ICrenshaw et al.l 120031 ; iKormendy fc Kennicutt 



2004 



Pierce et al.l 120071 : IShapiro et al.l 120081 ) . It is probable that if different forms of AGN ac- 
tivity are caused by contrasting fueling mechanisms, the relationship between various AGN 
spectral types may be different for the high and moderate luminosity populations. 

This study finds that the average type 1 and type 2 AGN hosts have similar stellar 
populations and similar levels of dust attenuation at z < 1. At z > 1 type 1 and type 2 AGN 
hosts appear to be fundamentally different. Even though the high redshift type 1 and type 2 
hosts have similar stellar populations, the type 1 AGN and host B band luminosity function 
requires that type 1 AGN hosts have significantly lower levels of dust extinction than what 
is required for type 2 AGN hosts to be in agreement with the observed near and mid-IR 
number counts. At z > 1, type 1 AGN hosts are intrinsically less dusty than type 2 AGN 
hosts. This suggests that the orientation based unified model works well for describing the 
local Seyfert population, but may not be appropriate for the high redshift, high luminosity 



quasar popu 
st udies (e.g ., 



Lutz et al 



ation. This conclusion is consistent with the find i ngs of several other recent 



Ballantvne et al.ll2006bl ; IHopkins fc Hernquistll2009l : iDraper &: Ballantyndl2010 



2010I ). For galaxies which are evolving secularly, the unified model appears to be 



an apt description. However, the violent growth experienced by black holes and their host 
bulges during major merger events does not appear to fit into the orientation based unified 
scheme. 

It is likely that different AGN fueling mechanisms will result in different relationships 
between AGN spectral types. In the merger scenario, young AGN ar e highly obscured 



and old quasars are unobscured (e.g., ISanders et al.l Il988l ; iFabianl Il999l : IPage et al.l I2004J : 



Hopkins et al.l 120061 ). Different spectral types of AGN fueled by secular processes related to 
nuclear starbursts are expected to be in agreement with the orientation based unified model 
due to the disky nature of th e nuclear starburst which is likely to both fuel and obscure the 
AGN (e.g., iBallantynd 120081 ). In the evolving model, it is found that mergers can play a 
strong role in fueling high Lx, z > 1 quasars. However, in the non-evolving model, when 
the enhanced star formation sources are considered, it appears that the dominant fueling 
mechanism is not mergers. Instead, the difference between spectral types may only be that 
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some galaxies have less dust and gas than others at z > 1. In the non-evolving far model, it 
appears that secular processes are the dominate AGN fueling mechanism at all redshift. For 
both the evolving and non-evolving model, at z < 1 the processes which lead to AGN activity 
are most likely secular. However, the findings of this study suggest that the dominant AGN 
fueling process changes at z ~ 1, since at z < 1 the average type 1 and type 2 AGN hosts are 
more similar than at higher redshift. In order to understand the mechanisms which trigger 
and fuel AGN, it is important for future studies to pay careful attention to which fueling 
mechanisms are dominant in different subsets of the AGN population. 

Thus, it has been shown that the relationship between different AGN spectral types 
is a helpful tool for understanding the dominant processes responsible for triggering and 
fueling AGN activity. It is found that at z < 1 the orientation based unified model holds, 
suggesting that the dominant AGN fueling mechanisms at z < 1 are secular processes. At 
z > 1, the orientation based unified model does not seem to hold. The evolving CT AGN 
model suggests that at z > 1 mergers are an important AGN fueling mechanism; however, 
the non-evolving CT AGN model suggests that even at z > 1 the dominant AGN fueling 
mechanisms are secular processes and that mergers play only a minor role in fueling AGN 
at all redshifts. 



4.5. Lx and redshift Evolution of f 2 

Understanding the evolution of ji is an important step in understanding the nature 
of AGN obscuration and the interplay between an AGN and it's host galaxy. If f 2 evolves 
with AGN L x , this su ggests that the o bscuring material is close enough to t he central engine 



that dust sublimation (JLawrencelll99ll ) and/or radiation pressure fed winds ( jKonigl fc Kartje 
19941 ) affect portions of the obscuring material. If f 2 evolves with redshift this would indicate 
that the evo lution of the obscuring m aterial is somehow connected to the evolution of the 
host galaxy (IBallantyne et al.ll2006bl ). While most model constraints used in this study are 



not very sensitive to the evolution of f 2 , the type 1 AGN and host B band luminosity function 
does prove to be a valuable test of the evolution of f 2 . This study finds that f 2 must evolve 
with AGN Lx and that f 2 probably evolves with redshift. 

If f 2 does not evolve with Lx, the predicted 1 < z < 2 type 1 AGN and host B band 
luminosity functions are not in agreement with observations. Using the host galaxy stellar 
population age and E(B — V) as free parameters does not allow for an appropriate fit to the 
data. The over-prediction of tub ~ 23 sources at 1 < z < 2 is not due to observational bias 
as rriB ~ 23 is significantly brighter than the depth of B band coverage ac cessible to surveys 



(e.g., B band 5a level for COSMOS is m B = 26.7; IScoville et al.l 120071 ). lLawrence &: Elvis 
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( 120 lOl ) argue that f 2 does not evolve with Lx and that apparent evolution of f 2 with Lx is 
due to X-ray observational biases. However, the z > 1 type 1 AG N and host B band luminos- 



ity function data u sed h ere is based on op t ical s election criteria (IHartwick fc Schaddll990 



Croom et al.ll2004j ). The iDella Ceca et al.l (119961 ) type 1 AGN and host B band luminosity 
function does use a sample defined through X-ray selection, and the z < 0.4 predicted type 1 
AGN and host B band luminosity function is in reasonable agreement with observation even 
if f 2 does not evolve with Lx- This suggests that X-ray observational biases do not create 
an artificial evolution of f 2 with Lx- Therefore, at least for z > 1, f 2 must evolve with Lx- 

While this study cannot rule out the possibility that f 2 is constant with redshift, our find- 
ings suggest that f 2 does e yolve mildly with redshift. While some studies su g gest that f 2 must 



evolve with redshift (e.g. lLa France et al.ll2005l ; iBallantyne et al.l l2006a| jbl; iTreister fc Urry 



20061 ). the evidence is still tentative. If AGN are fueled by different mechanisms, one might 
assume that the type 1/type 2 ratio will be different for the different mechanisms. Therefore, 
if different mechanisms dominate during different epochs, the type 1/type 2 ratio should also 
include some redshift evolution. It is possible that the slower secular evolution processes of 
Seyfert galaxies do not require evolution of f 2 with redshift or that the evolution is very 
mild, since this form of galaxy evolution is governed by stochastic processes. However, in 
the quasar regime, where AGN fueling is initiated by major mergers, a stronger evolution 
of f 2 with redshift may be necessary. Indeed, the type 1 AGN and host B band luminosity 
function which is least well fit by a non-evolving f 2 is the highest redshift bin considered, 
which is the redshift bin closest to the peak of quasar activity. While this study does not 
conclusively show that f 2 evolves with redshift, the findings presented here suggest that f 2 
does exhibit some evolution with redshift. 

It is found that f 2 evolves with Lx and f 2 is likely to evolve with redshift. This suggests 
that the obs curing medium is close enough to the central engine to be affected by dus t 
sublimation (JLawrencelll99ll ) and/or radiation pressure fueled winds (IKonigl fc Kartjdll994f ) . 
Also, the type 1 AGN and host B band luminosity function is an excellent X-ray independent 
tool to test the fraction of type 1 AGN predicted by X-ray observations. 



5. Summary 

By applying observational constraints at optical through mid-IR wavelengths to AGN 
host stellar population models, the average stellar properties of AGN hosts have been con- 
strained and hard X-ray through mid-IR SEDs have been developed which include emission 
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from the AGN and the host galaxy stellar population and dust enshrouded star formation.]. 
The findings of this study are summarized as follows. 

•Type 1 and type 2 AGN hosts have similar stellar populations at z < 1. 

•At z > 1, type 1 and type 2 AGN hosts have stellar populations of similar age, but 
type 1 AGN hosts are intrinsically less dusty than type 2 AGN hosts. 

•The orientation based unified model provides a good description of the z < 1 Seyfert 
population. The unified model does not seem to hold at z > 1, where quasar activity 
triggered by major mergers becomes more prevalent. 

•Multi- wavelength data is consistent with the paradigm in which (high Eddington ratio) 
CT AGN have recently undergone intense star formation. In the non-evolving model the 
AGN activity does not start until the starburst is >1 Gyr old and the stellar population 
is highly obscured by dust. In the evolving model, it is possible that the average high 
Eddington ratio CT AGN host galaxy is a LIRG or ULIRG, with the AGN activity starting 
one average >0.3 Gyr after the starburst. 

• >5-15% of the type 1 and type 2 AGN population may have enhanced levels of star 
formation, with average SFR « 20 M yr _1 . 

•If the enhanced star formation sources are included, the evolving model CT AGN 
hosts are unaffected, but the non-evolving model CT AGN hosts are not consistent with the 
AGN evolution scheme and instead are consistent with the orientation based unified scheme. 
Thus, if CT AGN are similar to type 2 AGN, then CT AGN host galaxies are similar to type 
2 AGN host galaxies; however, if CT AGN are a distinct population of AGN in a special 
evolutionary stage, then the host galaxies of CT AGN are also in a special evolutionary 
stage. 

•In order to fit the type 1 AGN and host B band luminosity function it is necessary 
for f 2 to evolve with Lx- However, future work studying the dependence of f 2 on redshift 
is necessary to elucidate the connection between AGN and their hosts. Understanding the 
redshift evolution of f% will offer insight into how AGN are fueled and the nature of AGN 
obscuration. 

•While the actual va lues of E(B — V) depend on the metallicity and IMF used in the 



Bruzual fc Charlotl (120031 ) host galaxy population synthesis model, the qualitative findings 
of this study are independent of, and robust against the uncertainties in, the host galaxy 
metallicity and IMF. 



2 AGN and host SEDs are available by contacting the authors. 
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Table 1. Ages of stellar populations in Gyr for various redshift bins. 

Type 1 and 2 AGN host Low Eddington ratio CT AGN host 1 

YSP f DSP* YSP+ DSP* 

(0.1M*) (0.9M*) (0.1M*) (0.9M*) 



z < 0.5 


2.0 


7.0 


4.5 


10 


0.5 < z < 1 


2.0 


7.0 


4.5 


7.0 


1 < ^ < 2 


2.0 


4.5 


4.5 


4.5 


2 < z < 3 


2.0 


2.0 


2.0 


2.0 


3 < z 


1.0 


1.0 


1.0 


1.0 



^The Low Eddington ratio CT AGN hosts are only relevant in the evolving model. The 
age of the stellar population of the evolving model high Eddington ratio CT AGN hosts is 
the same in all z bins and is shown in the bottom row of Table [3j 



^Younger Stellar Population 
^Dominant Stellar Population 



Table 2. E(B — V) for the different IMFs and metallicities considered. The first row 
summarizes the working model of this paper. 

Model Unified/type 2 AGN host Type 1 AGN host 

IMF Metallicity z < 1 z > 1 z<\ z > 1 



Chabrier (20031 


Z & 


0.25 


0.5 


0.25 


0.05 


Chabrier ('2003) 


0.2Z Q 


0.3 


0.5 


0.3 


0.1 


Chabrier (2003) 


2.5Z Q 


0.2 


0.7 


0.1 


0.05 


SalDeter (1955) 


Z & 


0.1 


0.1 


0.1 


0.0 
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Table 3. Summary of the AGN host stellar population parameters for the working model 
of this study. The SFRs in parenthesis refer to the average SFR of the enhanced star 

formation sources. 



AGN population 



Stellar population age E(B-V) E(B-V) SFR 

(Gyr) z < 1 z > 1 (M yr" 1 ) 



Type 1 AGN 
Unified/type 2 AGN 



col 2 & 3 Table [TJ 0.25 0.05 2.0(20) 

col 2 & 3 Table [TJ 0.25 0.5 2.0(20) 



Non-evolving model CT AGN 

w/o enhanced star formation sources 
w/ enhanced star formation sources 



1.0 1.0 

col 2 & 3 Table [2 0.25 



1.0 


10 


0.5 


2.0 (20) 


1.0 


1.0 


0.6 


175 



Evolving model CT AGN 

Low Eddington ratio CT AGN 
High Eddington ratio CT AGN 



col 4 & 5 Table [2 1.0 

0.3 0.6 
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Fig. 1. — Smoothed AGN and host rest frame SEDs. Solid lines sh ow the total AGN 



and h ost SED. Dot-dashed lines show AGN SEDs, dashed lines show iBruzual fc Chariot 
(120031 ) stellar population model SEDs, and dotted lines show iRieke et al.l ( 120091 ) dusty star 
formation templates, (a) type 1 AGN with Lx = 10 43 erg s _1 , a stellar population of M* 
= 10 11 M with 90% of M* in a 4.5 Gyr old stellar population and 10% of M* in a 2 Gyr 
old stellar population with E(B — V) « 0.25, and the L/# = 10 10 L dusty star formation 
template, corresponding to SFR s=s 2 M yr _1 . (b) type 2 AGN with Lx = 10 43 erg s _1 , 
the same stellar population as in (a) but with E(B — V) ~ 0.50, and the Lir = 10 11 L Q 
template, corresponding to SFR w 17 M yr _1 . (c) low Eddington ratio CT AGN with Lx 
= 10 42 erg s" 1 , a 4.5 Gyr old stellar population of M* = 10 12 M Q with E(5 - V) « 1.0, 
and SFR w 1 M Q yr" 1 , corresponding to the L IR = 10 975 L template, (d) high Eddington 
ratio CT AGN with L x = 10 44 erg s -1 , a 0.3 Gyr old stellar population of M* = 10 10 M 
with E(B - V) « 0.50, and the L JB = 10 12 L template, or SFR w 175 M yr" 1 . 
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Fig. 2. — J band space density for AGN and hosts for the unified AGN host model. The black 
lines show the AGN and host J band space density for the model host galaxies described in 
row 2 of Table El The solid lines show the total AGN and host J band space density while 
the dotted lines show the space density for type 1 AGN and hosts, dashed lines show the 
space density for type 2 AGN and hosts, and dot-dashed lines show th e CT AGN and h ost 
space density. Data from the mid-IR and X-ray selected AGN sample of lAssef et al.l ( 120111 ) is 
also shown. The obvious discrepancy at z < 0.75 between the model presented here and the 
ob servations reported by lAssef et al.l (1201 ll ) are primarily due to a selection criterion used 
by lAssef et al.l (120 111 ) for optically extended sources at z < 0.75 which cannot be replicated 
here as the models do not contain information of the spatial extent of the host galaxies. 
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Fig. 3. — Soft X-ray flux versus R band magnitude for AGN and hosts for the non-evolving 
model. AGN and hosts are shown for L x < 10 45 erg s _1 , z < 3, and M* = 10 95 , 10 10 , 
10 10 - 5 , 10 11 , 10 1L5 , and 10 12 M . The green triangles show the unified AGN hosts model 
and type 2 AGN and hosts, blue circles show type 1 AGN and hosts, and red squares show 
non-evolving model CT AGN and hosts. The model host galaxies shown here are described 
by the parameters in the first three rows of Table [3J The dashed lines show the empirical 
relationship log(/ .5-2fcey//R) — 0±1. The horizontal dotted line marks R = 25.5, above 
which the source is considered an optically faint X-ray AGN. 
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Fig. 4. — Type 1 AGN and host B band luminosity function. Different colors show the 
luminosity function in different redshift bins. Black lines and data refer to z < 0.4, blue 
lines and data refer to 1.0 < z < 1.55, and red lines and data refer to 1.55 < z < 2.1. The 
blue and red dashed lines show the luminosity function for the unified AGN host model as 
summarized in row one of Table [3j The black dashed line shows the dustiest average AGN 
host at z < 1 allowed by the unified AGN host model, E(B — V) ~ 0.4. The solid lines 
show the type 1 AGN host best fit model as summarized in row 2 of Table [3j The AGN 
contribution to the luminosity function is shown as the dot-dashed lines. Data points show 
vario us type 1 AGN and hos t B band luminosity functions from the literature: squares are 
from iDella Ceca et al.l (Il996f ). triangles are from ICroom et all (J2004f ). and stars are the qo 
= 0.5 luminosity function from lHartwick fc Schadd ( )1990l ) converted to the cosmology used 
here. 
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Fig. 5. — Near and mid-IR number counts for AGN and hosts for the non-evolving model 
with /2_8 > 1 x 10~ 16 erg s" 1 cm -2 . The cyan lines show the predicted number counts 
for AGN alone. The black lines show the AGN and host number counts for the model host 
galaxies described by the first three rows of Table [3j The solid lines show the total AGN and 
host number counts while the dotted lines show the counts for type 1 AGN and hosts, dashed 
lines show the counts for type 2 AGN and hosts, and dot-dashed lines show the CT AGN 
and h ost counts. Data is from Spitzer observations of X-ray selected AGN in the GOODS 
fields JTreister et all bpOfih . 
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Fig. 6. — J band space density for AGN and hosts for the non-evolving model. The black 
lines show the AGN and host J band space density for the model host galaxies described in 
the first three rows of Table O The solid lines show the total AGN and host J band space 
density while the dotted lines show the space density for type 1 AGN and hosts, dashed lines 
show the space density for type 2 AGN and hosts, and dot-dashed lines show the CT AGN 
and host space density. Data points are the same as in Figure [2j 
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Fig. 7. — Near and mid-IR number counts for AGN and hosts for the evolving model with 
/2_8 > 1 x 10~ 16 erg s _1 cm -2 . As in Figure the cyan solid lines show the predicted 
number counts for AGN alone, the black solid lines show the total AGN and host number 
counts, the black dotted lines show the counts for type 1 AGN and hosts, and the black 
dashed lines show the counts for type 2 AGN and hosts. The type 1 and type 2 AGN hosts 
are the same as in Figure |5j The black dot-dashed lines show the total CT AGN and host 
number counts for the evolving model. The low Eddington ratio CT AGN hosts are shown 
by the blue dot-dashed lines and the high Eddington ratio CT AGN hosts are shown by the 
red dot-dashed lines. The evolving model CT AGN host galaxies are summarized in Table 
[3j The data shown is the same as in Figure [51 
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Fig. 8. — J band space density for AGN and hosts for the evolving model. The solid lines 
show the total AGN and host J band space density. The dotted lines show the type 1 AGN 
and hosts space density and the dashed lines show the type 2 AGN and hosts space density. 
The type 1 and type 2 AGN hosts have the same stellar populations as in Figure [6j The black 
dot-dashed lines show the total CT AGN and host J band space density for the evolving 
model. The low Eddington ratio CT AGN hosts, shown by the blue dot-dashed lines, and 
the high Eddington ratio CT AGN hosts, shown by the red dot-dashed lines, are described 
in Table |3j Data points are the same as in Figure [2j 
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Fig. 9. — Soft X-ray flux versus R band magnitude for AGN and hosts for the evolving 
model. AGN and hosts are shown for Lx < 10 45 erg s _1 , z < 3, and M* = 10 9 ' 5 , 10 10 , 
10 10 ' 5 , 10 11 , 10 1L5 , and 10 12 M . The blue circles show type 1 AGN and hosts and green 
triangles show type 2 AGN and hosts. The type 1 and type 2 AGN hosts are the same as 
in Figure [3j The red squares show CT AGN and hosts. The red filled squares show the low 
Eddington ratio CT AGN and hosts, while the open red squares show the high Eddington 
ratio CT AGN and hosts. The evolving model CT AGN host galaxies are described in Table 
El The horizontal dotted line again marks R = 25.5, above which the source is considered 
an optically faint X-ray AGN. 
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Fig. 10. — Type 1 AGN and host B band luminosity function when f2 is assumed to not 
evolve with z. Colors and data points are the same as in Figure HI The z < 1 type 1 AGN 
hosts have E(B — V) w 0.25 and at higher redshift the type 1 AGN hosts have E(B — V) « 
0.4, in contrast to the case where f% does evolve with z and type 1 AGN hosts have a lower 
E[B — V) at higher z than locally. 
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Fig. 11.— i? versus f 24 for evolving model AGN with L x = 10 42 , 10 43 - 5 , 10 45 and 10 465 erg 
s _1 , z < 5, and M* = 10 10 , 10 11 , and 10 12 M for various soft X-ray flux ranges. Panel a 
shows all AGN regardless of fo.5-2, panel b shows all AGN with /0.5-2 < 10"~ 15 erg s _1 cm" 2 , 
and panel c shows all AGN with fo.5-2 < 10~ 16 erg s _1 cm" 2 . The blue circles show the type 
1 AGN and the green triangles show the type 2 AGN. The open points mark the enhanced 
star formation objects. The red filled squares show the low Eddington ratio CT AGN and 
the open red squares show the high Eddington ratio CT AGN. The thick black line shows 
where /24//H = 1000, therefore the area of interest is above the thick line. The vertical lines 
mark / 24 = 550 /iJy and / 2 4 = 700 /xJy. 
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Fig. 12.— Redshift distribution of R versus / 24 for CT AGN with L x = 10 42 , 10 43 - 5 , 10 45 
and 10 46 ' 5 erg s _1 and M* = 10 10 , 10 11 , and 10 12 M using the evolving model for various 
soft X-ray flux ranges. Panel a shows all CT AGN regardless of /0.5-2, panel b shows the CT 
AGN with /o. 5 _2 < 10~ 15 erg s" 1 cm" 2 , and panel c shows the CT AGN with / . 5 _ 2 < 10~ 16 
erg s _1 cm" 2 . Point styles designate AGN in different redshift ranges. Blue circles show 
AGN with z < 1, green triangles show AGN with 1 < z < 2, red squares show AGN with 2 
< z < 3, and cyan stars show AGN with z > 3. The filled points show the low Eddington 
ratio CT AGN and the open points show the high Eddington ratio CT AGN. The thick black 
line shows where fi±/ Jr — 1000, therefore the area of interest is above the thick line. The 
vertical lines mark f 24 = 550 /xJy and f 24 = 700 /zJy. 



